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using a time-dependent Dalitz plot analysis. The results are obtained from a data sample of 347 
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B Factory at SLAC. We measure 26 coefficients of the bilinear form factor terms occurring in the 
I time-dependent decay rate of the B^ meson and derive the physically relevant quantities from these 
coefficients. In particular we find a three standard deviation evidence of direct CP- violation in the 
B^ —I- p^TT^ decays, with systematic uncertainties included. We also achieve a constraint of the 
angle a of the Unitarity Triangle. All results presented are preliminary. 
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1 INTRODUCTION 



Measurements of the parameter sin2/3 [1, 2] have estabUshed CP violation in the meson system 
and provide strong support for the Kobayashi and Maskawa model of this phenomenon as arising 
from a single phase in the three-generation CKM quark-mixing matrix [3]. We present, in this paper, 
results from a time-dependent analysis of the — >■ tt+tt^tt'^ Dalitz plot (DP) that is dominated by 
the p{770) intermediate resonances of all charges and their interference. The goal of the analysis 
is the simultaneous extraction of the strong transition amplitudes and the weak interaction phase 
a = arg [—ytdYtb/^udYubl Unitarity Triangle. In the Standard Model, a non-zero value for 

a is responsible for the occurrence of mixing-induced CP violation in this decay. The BABAR and 
Belle experiments have obtained constraints on a from the measurement of effective quantities 
sin2aefr in B decays to tt+tt" [4, 5] and from B decays to p'^p~ [6, 7], using an isospin analysis [8]. 

Unlike n^n~, p^ir^ is not a CP eigenstate, and four flavor-charge configurations {B^{B^) — > 
p^TT^) must be considered. The corresponding isospin analysis [9] is unfruitful with the present 
statistics since two pentagonal amplitude relations with 12 unknowns have to be solved (compared 
to 6 unknowns for the vr+vr^ and p^p" systems). However, it has been pointed out by Snyder 
and Quinn [10], that one can obtain the necessary degrees of freedom to constrain a without 
ambiguity by explicitly including in the analysis the variation of the strong phases of the interfering 
p resonances in the Dalitz plot. 



1.1 DECAY AMPLITUDES 

We consider the decay of a spin-zero B^ with four-momentum p b into the three daughters vr"*" , vr" , 
7r°, with and po their corresponding four- momenta. Using as independent (Mandelstam) 

variables the invariant squared masses 

s+ = ip++Pof , s- = ip-+pof , (1) 

the invariant squared mass of the positive and negative pion, sq = {p+ +P-)'^, is obtained from 
energy and momentum conservation 

So = rn%o + 2m^+ -|- m^o - s+ - s_ . (2) 

The differential B^ decay width with respect to the variables defined in Eq. (1) (i.e., the Dalitz 
plot) reads 

<ir(B"^.W)= pl3^&,&_, (3) 

where ^3,^ is the Lorentz-invariant amplitude of the three-body decay. 

We assume in the following that the amplitudes A^,^ and its complex conjugate As-j^, corre- 
sponding to the transitions B° — Tr+Tr^Tr*^ and B^ — > tt+tt^tt^, respectively, are dominated by the 
three resonances p+, p- and /9°. The p resonances are assumed to be the sum of the ground state 
/j(770) and the radial excitations p(1450) and /9(1700), with resonance parameters determined by 
a combined fit to r"*" UtTt'^'k^ and e^e~ — > 7r"'"7r~ data [11]. Since the hadronic environment is 
different in B decays, we cannot rely on this result and therefore determine the relative p(1450) and 
p(1700) amplitudes simultaneously with the CP parameters from the fit. Variations of the other 
parameters and possible contributions to the B^ — ^ 7r"'"7r~7r'^ decay other than the p's are studied 
as part of the systematic uncertainties (Section 4) . 
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Including the B^B^ mixing parameter q/p into the B^ decay amphtudes, we can write [10, 12] 

As. = + f_A- + foA\ (4) 

As^ = UA+ + f_A- + foA\ (5) 

where the (with k = {+, —,0} denote the charge of the p from the decay of the S° meson) 
are functions of the DaHtz variables s+ and s_ that incorporate the kinematic and dynamical 
properties of the B^ decay into a (vector) p resonance and a (pseudoscalar) pion, and where the 
A'^ are complex amplitudes that may comprise weak and strong transition phases and that are 
independent of the Dalitz variables. Note that the definitions (4) and (5) imply the assumption 
that the relative phases between the p(770) and its radial excitations are CP-conserving. 

Following Rcf. [11], the p resonances are parameterized in by a modified rclativistic Breit- 
Wigner function introduced by Gounaris and Sakurai (GS) [15]. Due to angular momentum con- 
servation, the spin-one p resonance is polarized in a helicity-zero state. For a p'^ resonance with 
charge k, the GS function is multiplied by the kinematic function — 4|p^||pt| cos ^k, where is 
the momentum of either of the daughters of p-resonancc defined in the p-resonance rest frame, and 
where is the momentum of the particle not from p decay defined in the same frame, and cos 6^ 
the cosine of the helicity angle of the p'^. For the p"*" (p~), 9+ {9-) is defined by the angle between 
the vr*' (tt^) in the p"*" (p~) rest frame and the p^ (p~) flight direction in the B^ rest frame. For 
the p^, 9o is defined by the angle between the vr"'" in the p^ rest frame and the p^ flight direction 
in the B^ rest frame. With these definitions, each pair of GS functions interferes destructively at 
equal masses-squared. 

The occurrence of cos 0^ in the kinematic functions substantially enhances the interference be- 
tween the different p bands in the Dalitz plot, and thus increases the sensitivity of this analysis [10]. 

1.2 TIME DEPENDENCE 

With At = tsTf — ttag defined as the proper time interval between the decay of the fully reconstructed 
and that of the other meson B^^^^, the time-dependent decay rate |^^^(At)p (|^3^(At)|^) when 
the tagging meson is a B^ (B°) is given by 

|^37rP+|^37rp¥(|^3^P - l^aTrP) cos(AmdAt)± 2Im [^3^A^^] sin(Am<iAt) 

(6) ■ 

where r^o is the mean B^ lifetime and Arrid is the B^B^ oscillation frequency. Here, we have 
assumed that CP violation in B^B^ mixing is absent {\q/p\ = 1), AFg^ = and CPT is conserved. 
Inserting the amplitudes (4) and (5) , one obtains for the terms in Eq. (6) 

kG{+,-,0} K<crG{ + ,-,0} 

Im(^3.Ay = + E (Re[A/*]/^+Im[/,/;]e) , (7) 

Ke{+,-,o} K;<<Te{+,-,o} 
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with 



ut = Ivl'^P ± lA'^P , (8) 
jj±,Ke(Im) ^ Rc(Im) [A^'A"* ± A" A"*] , (9) 

= lm[A''A''*] , (10) 

I^^ = Re [I'^A'^* - A^A'^*] , (11) 

7;^™ = Im + A'^^'^*] . (12) 

The 27 coefficients (8)-(12) are real-valued parameters that multiply the /«/* bilinears (where 

K and cj denote the charge of the p resonances) [16]. These are the observables that arc determined 
by the fit. Each of the coefficients is related in a unique way to physically more intuitive quantities, 
such as tree-level and penguin-type amplitudes, the angle a, or the quasi-two-body CP and dilution 
parameters [17] (c/. Section 6). The parameterization (7) is general; the information on the mirror 
solutions (e.g., on the angle a) that are present in the transition amplitudes A'^, A'^ is conserved. 
In this paper, we determine the relative values of U and I coefficients to U^. 

The choice to fit for the U and I coefficients rather than fitting for the complex transition 
amplitudes and the weak phase a directly is motivated by the following technical simplifications: 
(i) in contrast to the amplitudes, there is a unique solution for the U and / coefficients requiring only 
a single fit to the selected data sample, (ii) in the presence of background, the U and / coefficients 
are approximately Gaussian distributed, which in general is not the case for the amplitudes, and 
{in) the propagation of systematic uncertainties and the averaging between different measurements 
are straightforward for the ?7's and /'s. 

The coefficients are related to resonance fractions (branching fractions and charge asym- 
metries); the determine the relative abundance of the decay into p'^7r~ and p~7r~^ and the 
time-dependent direct CP asymmetries. The 7^ measure mixing-induced CP violation and are sen- 
sitive to strong phase shifts. Finally, the C/^^'^^^™'* and 7^''^™^ describe the interference pattern 
in the Dalitz plot, and their presence distinguishes this analysis from the quasi-two-body analysis 
previously reported in [17]. They represent the additional degrees of freedom that allow one to 
determine the unknown penguin pollution and the relative strong phases. However, because the 
overlap regions of the resonances are small and because the events reconstructed in these regions 
suffer from large misreconstruction rates and background, a substantial data sample is needed to 
perform a fit that constrains all amplitude parameters. 

We determine the quantities of interest in a subsequent least-squares fit to the measured U and 
7 coefficients. 

1.3 NORMALIZATION 

The decay rate (6) is used as a probability density function (PDF) in a maximum-likelihood fit and 
must therefore be normalized: 

\AUAt)\' I ,J AUm' , (13) 



where 



(1^3.1' + 1^3.!') = E m')U: + 2ReJ2 {f.K){Ut/'' + i-U+/^) . (14) 

Ke{+,-,Q} K<ae{+,-,Q} 
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The complex expectation values {f^fa) obtained from high-statistics Monte Carlo integration 
of the Dalitz plot (3), taking into account acceptance and resolution effects. 

The normalization of the decay rate (6) renders the normalization of the U and I coefficients 
arbitrary, so that we can fix one coefficient. By convention, we set = 1. 

1.4 THE SQUARE DALITZ PLOT 

Both the signal events and the combinatorial e+e~ — > qq {q = u,d,s, c) continuum background 
events populate the kinematic boundaries of the Dalitz plot due to the low final state masses 
compared to the mass. We find the representation (3) is inadequate when one wants to use 
empirical reference shapes in a maximum- likelihood fit. Large variations occurring in small areas 
of the Dalitz plot are very difficult to describe in detail. These regions are particTilarly important 
since this is where the interference, and hence our ability to determine the strong phases, occurs. 
We therefore apply the transformation 

ds+ ds- — > I det J| dm' dd' , (15) 

which defines the Square Dalitz plot (SDP). The new coordinates arc 



1 / niQ — m™™ \ _ 1 



2 . - 1 , ^' ^ -^0, (16) 



m = — arccos 

TT \ UIq"^ — mfj"" J IT 

where mo is the invariant mass between the charged tracks, m^^^ = m^o — m^o and rn™"^ = 2m^+ 
are the kinematic limits of mo, ^o is the helicity angle, and J is the Jacobian of the transformation 
that zooms into the kinematic boundaries of the Dalitz plot. The new variables range between 
and 1. The determinant of the Jacobian is given by 

I , , ^1*11*1 ^"^0 5cos6'o , . 

IdetJl = 4|p+||Po|mo-^-^^, (17) 

where \p*^_\ = — and |po| = Eq — m?o , and where the energies E^ and Eg are in 

the TT+TT^ rest frame. Figure 1 shows the determinant of the Jacobian as a function of the SDP 
parameters m' and 9'. If the events in the nominal Dalitz plot were distributed according to a 
uniform (non-resonant) prior, their distribution in the SDP would match the plot of | det J\. 

The effect of the transformation (15) can be seen by looking at Figure 2 which displays the 
nominal and square Dalitz plots for signal events generated with toy Monte Carlo (MC): the 
homogenization of the distribution is clearly visible. This simulation does not take into account any 
detector effect and corresponds to a particular choice of the decay amplitudes for which destructive 
interferences occur at equal p masses. To simplify the comparison, hatched areas showing the 
interference regions between p bands and dashed isocontours y/sZjTZ^ = 1.5 GeV/c^ have been 
superimposed on both Dalitz plots. 



2 THE BABAR DETECTOR AND DATASET 

The data used in this analysis were collected with the BABAR detector at the PEP-II asymmetric- 
energy e~^e~ storage ring at SLAC between October 1999 and June 2006. The sample consists of 
about 310 fb~\ corresponding to (346 ± 3) x 10^ BB pairs collected at the T(45') resonance ("on- 
resonance"), and an integrated luminosity of 21.6 fb~^ collected about 40 MeV below the T{AS) 
( "off-resonance" ) . 



11 



A detailed description of the BABAR detector is presented in Ref. [18]. The tracking system used 
for track and vertex reconstruction has two main components: a sihcon vertex tracker (SVT) and a 
drift chamber (DCH), both operating within a 1.5 T magnetic field generated by a superconducting 
solenoidal magnet. Photons are identified in an electromagnetic calorimeter (EMC) surrounding 
a detector of internally reflected Chercnkov light (DIRC), which associates Cherenkov photons 
with tracks for particle identification (PID). Muon candidates are identified with the use of the 
instrumented flux return (IFR) of the solenoid. 

3 ANALYSIS METHOD 

The U and / coefficients and the Tr+Tr^Tr" event yield are determined by a maximum- 

likelihood fit of the signal model to the selected candidate events. Kinematic and event shape 
variables exploiting the characteristic properties of the events are used in the fit to discriminate 
signal from background. 

3.1 EVENT SELECTION AND BACKGROUND SUPPRESSION 

We reconstruct — > 7r~^7T~7T^ candidates from pairs of oppositely-charged tracks, which are re- 
quired to form a good quality vertex, and a 7r° candidate. In order to ensure that all events are 
within the Dalitz plot bounaries, we constrain the three-pion invariant mass to the B-mass. We 
use information from the tracking system, EMC, and DIRC to remove tracks for which the PID 
is consistent with the electron, kaon, or proton hypotheses. In addition, we require that at least 
one track has a signature in the IFR that is inconsistent with the muon hypothesis. The tt*^ can- 
didate mass must satisfy 0.11 < 771(77) < 0.16GeV/c^, where each photon is required to have an 
energy greater than 50 McV in the laboratory frame (LAB) and to exhibit a lateral profile of energy 
deposition in the EMC consistent with an electromagnetic shower. 

A S-meson candidate is characterized kinematically by the energy-substituted mass mss = 
[(^s -I- po • Pb)'^/Eq — p^]2 and energy difference AE = — \-y/s, where {Eb,'Pb) and (£^0)Po) 
are the four-vectors of the S-candidate and the initial electron-positron system, respectively. The 
asterisk denotes the T(46') frame, and s is the square of the invariant mass of the electron-positron 
system. We require 5.272 < mES < 5.288 GeV/c^, which retains 81% of the signal and 8% of the 
continuum background events. The AE resolution exhibits a dependence on the tt^ energy and 
therefore varies across the Dalitz plot. We account for this effect by introducing the transformed 
quantity AE' = {2AE-AE+-AE_)/{AE+-AE_), with AE±{mo) = c±-{c±^c) {mo/m^''''y , 
where mo is strongly correlated with the energy of 7r°. We use the values c = 0.045 GeV, c_ = 
-0.140 GeV, c+ = 0.080 GcV, = 5.0 GeV, and require -1 < AE' < 1. These values have 

been obtained from Monte Carlo simulation and are tuned to maximize the selection of correctly 
reconstructed over misreconstructed signal events. The requirement retains 75% (25%) of the signal 
(continuum). 

Backgrounds arise primarily from random combinations in continuum events. To enhance dis- 
crimination between signal and continuum, we use a neural network (NN) [19] to combine four 
discriminating variables: the angles with respect to the beam axis of the B momentum and B 
thrust axis in the T{AS) frame, and the zeroth and second order monomials Lo,2 of the energy flow 
about the B thrust axis. The monomials arc defined by Lj = X^jPi x Icos^jp, where 9i is the 
angle with respect to the B thrust axis of track or neutral cluster i, pi is its momentum, and the 
sum excludes the B candidate. The NN is trained in the signal region with off-resonance data and 
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simulated signal events. The final sample of signal candidates is selected with a requirement on 
the NN output that retains 77% (8%) of the signal (continuum). 

The time difference At is obtained from the measured distance between the z positions (along 
the beam direction) of the B^^ and B^^^ decay vertices, and the boost P'-f = 0.56 of the e~^e~ 
system: At = Az/(3^c. To determine the flavor of the B^^^ we use the B flavor tagging algorithm 
of Ref. [22]. This produces six mutually exclusive tagging categories. We also retain untagged 
events in a seventh category to improve the efficiency of the signal selection and because these 
events contribute to the measurement of direct CP violation. Events with multiple B candidates 
passing the full selection occur in 16% (p'^iT^) and 9% (p'^vr^) of the cases. If the multiple candidates 
have different tt^'s, we choose the candidate with the reconstructed 7r° mass closest to the nominal 
one; in the case that both candidates have the same 7r°, we pick the first one. 

The signal efficiency determined from MC simulation is 24% for B^ — > p^ir^ and B^ — p^ir^ 
events, and 11% for non-resonant B^ — > tt'^tt^'k^ events. 

Of the selected signal events, 22% of B° p^tt^, 13% of B^ P°7r°, and 6% of non-resonant 
events are misreconstructed. Misreconstruced events occur when a track or neutral cluster from 
the tagging B is assigned to the reconstructed signal candidate. This occurs most often for low- 
momentum tracks and photons and hence the misreconstructed events concentrate in the corners of 
the Dalitz plot. Since these are also the areas where the p-mesons overlap strongly, it is important 
to model the misreconstruced events correctly. The details of the model for misreconstructed events 
over the Dalitz plot is detailed in Section 3.3.1. 

3.2 BACKGROUND FROM OTHER B DECAYS 

We use MC simulated events to study the background from other B decays. More than a hundred 
channels have been considered in the preliminary studies, of which twenty-nine have been finally 
included in the likelihood model - decays with at least two events expected after selection. These 
exclusive i?-background modes are grouped into eighteen different classes gathering decays with 
similar kinematic and topological properties: six for charged charmless i?-decays, eight for neu- 
tral charmless i?-decays and four for exclusive neutral charmed i?-decays. Two additional classes 
account for inclusive neutral and charged 6 — > c decays. 

Table 1 summarizes the twenty background classes which are used in the fit. For each mode, 
the expected number of selected events is computed by multiplying the selection efficiency (esti- 
mated using MC simulated decays) by the branching fraction scaled up to the dataset luminosity 
(310 fb-^). The world average branching ratios have been used for the experimentally known de- 
cay modes. When only upper limits are given, they have been translated into branching ratios 
using all information available such as additional conservative hypotheses (e.g. 100% longitudinal 
polarization for B ^ pp decay) if needed. 

3.3 THE MAXIMUM-LIKELIHOOD FIT 

We perform an unbinned extended maximum-likelihood fit to extract the inclusive B^ — ^ 7r"^7r~7r^ 
event yield and the U and / coefficients defined in Eqs. (8)— (12). The fit uses the variables 
rTT-ESi AE' , the NN output, and the Dalitz plot to discriminate signal from background. The At 
measurement allows to the determination of mixing-induced CP violation and provides additional 
continuum-background rejection. 

The selected on-resonance data sample is assumed to consist of signal, continuum-background 
and B-background components, separated by the fiavor and tagging category of the tag side B 
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decay. The signal likelihood consists of the sum of a correctly reconstructed ( "truth- matched" , 
TM) component and a misreconstructed ("self-cross-feed", SCF) component. 

The probability density function (PDF) for an event i in tagging category c is the sum of 
the probability densities of all components, namely 



Vt = -^Stt/stt (1 - /sCF)'P37r-TM,i + /sCf'^Stt-SCF,' 

class -j^ 



Nf 

class 



+ E ^BOjrBOjno,ij , (18) 

where: N^-^ is the total number of 7r"'"7r~7r^ signal events in the data sample; is the fraction 
of signal events that are tagged in category c; /scf is the fraction of SCF events in tagging 
category c, averaged over the Dalitz plot; ^|^_tm i ^Stt-scf i products of PDFs of the 

discriminating variables used in tagging category c for TM and SCF events, respectively; A'^^ is 
the number of continuum events that are tagged in category c; gtag,i is the tag flavor of the event, 
defined to be for a B^^^^ and —1 for a B^^^, A^^^tag parameterizes possible tag asymmetry in 
continuum events; Vg^ i is the continuum PDF for tagging category c; N^^^ i^Sass) is number 
of charged (neutral) S-rclatcd background classes considered in the fit; N^+j (N^oj) is the number 
of expected events in the charged (neutral) S-background class j; f%+j i fsOj) is the fraction of 
charged (neutral) ^-background events of class j that are tagged in category c; ^B+,tagj describes 
a possible tag asymmetry in the charged-B background class j; correlations between the tag and 
the position in the Dalitz plot (the "charge") are absorbed in tag-flavor-dependent Dalitz plot 
PDFs that are used for charged-S and continuum background; is the i3^-background PDF 

for tagging category c and class j; finally, V%o is the neutral- i?-background PDF for tagging 
category c and class j. 

The PDFs {X = {TM, SCF, continuum, B - bkg}) are the product of the four PDFs of the 
discriminating variables, xi = mES, X2 = AE', xs = NNoutput, and the triplet X4 = {m', 9', At}: 

4 

- U^xmi^^)- (19) 

k=l 

The extended likelihood over all tagging categories is given by 

7 _ N" 

c=l i 

where A^'^ is the total number of events expected in category c. 

A total of 68 parameters, including the inclusive signal yield and the parameters from Eq. (6), 
are varied in the fit. Most of the parameters describing the continuum distributions are also floated 
in the fit. 



14 



3.3.1 THE At AND DALITZ PLOT PDFS 



The Dalitz plot PDFs require as input the Dalitz plot-dependent relative selection efficiency, e = 
€{m',6'), and SCF fraction, /scF = fscF{n^' ,9')- Both quantities are taken from MC simulation. 
Away from the Dalitz plot corners the efficiency is uniform, while it decreases when approaching the 
corners, where one out of the three bodies in the final state is close to rest so that the acceptance 
requirements on the particle reconstruction become restrictive. Combinatorial backgrounds and 
hence SCF fractions are large in the corners of the Dalitz plot due to the presence of soft neutral 
clusters and tracks. 

For an event i, we define the time-dependent Dalitz plot PDFs 

P3^-TM,i = Siil- fscF,i) I det Ji\ \Ai^{At)\^ , (21) 
P3n-SCF,i = eifscF,i\detJi\\Ai,{At)W (22) 

where PsTr-TM,^ and PsTr-SCF.i arc normalized. The corresponding phase space integration involves 
the expectation values (e(l'- /scf) I det J| /'^/'^*) and (e /scF I det J| /'^/'^*) for TM and SCF 
events, where the indices k, a run over all resonances belonging to the signal model. The expectation 
values are model-dependent and are computed with the use of MC integration over the square Dalitz 
plot: 

, . r„Vn £(l-/sCF)|det J|/«r*dmW , , 

£ 1 - /scF det J = -^0 Jo I 23 

^ ^ ^' ^ ' e \ det J\f-f-* dm' d9' ^ ' 

and similarly for (e | det J| f^f^*), where all quantities in the integrands are Dalitz plot-dependent. 

Equation (18) invokes the phase space-averaged SCF fraction /scf = (/scF | det J\ f^f"*). The 
PDF normalization is decay-dynamics-dependent and is computed iteratively. We determine the 
average SCF fractions separately for each tagging category from MC simulation. 

The width of the dominant /9(770) resonance is large compared to the mass resolution for TM 
events (about 8MeV/c^ core Gaussian resolution). We therefore neglect resolution effects in the 
TM model. Misreconstructed events have a poor mass resolution that strongly varies across the 
Dalitz plot. It is described in the fit by a 2 x 2-dimensional resolution function 

RscF{m',X,m[A) , (24) 

which represents the probability to reconstruct at the coordinate {m'^,6'j.) an event that has the 
true coordinate {m^,9[). It obeys the unitarity condition 

1 1 

j J RQcv{rnj.,6'j.,m[,0[) dm^dO'^ = 1, (25) 



and is convolved with the signal model. The -RsCF function is obtained from MC simulation. 

We use the signal model described in Section 1.1. It contains the dynamical information and 
is connected with At via the matrix element (6), which serves as PDF. It is diluted by the effects 
of mistagging and the limited vertex resolution [17]. The At resolution function for signal and 
5-background events is a sum of three Gaussian distributions, with parameters determined by a 
fit to fully reconstructed decays [22]. 

The Dalitz plot- and At-dependent PDFs factorize for the charged-S-background modes, but not 
(necessarily) for the neutral-B background due to B^B^ mixing. 
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The charged Z?-background contribution to the hkeUhood (18) involves the parameter ^s+^tag) 
multiphed by the tag flavor gtag of the event. In the presence of significant tag- "charge" correla- 
tion (represented by an effective flavor-tag- versus-Dalitz-coordinate correlation), it parameterizes 
possible direct CP violation in these events. Wc also use distinct square Dalitz plot PDFs for each 
reconstructed B flavor tag, and a flavor-tag-averaged PDF for untagged events. The PDFs are 
obtained from MC simulation and are described with the use of non-parametric functions. The 
At resolution parameters are determined by a flt to fully reconstructed S+ decays. For each S+- 
background class we adjust effective lifetimes to account for the misreconstruction of the event that 
modifies the nominal At resolution function. 

The neutral-S background is parameterized with PDFs that depend on the flavor tag of the 
event. In the case of CP eigenstates, correlations between the flavor tag and the Dalitz coordinate 
are expected to be small. However, non-CP eigenstates, such as af-K^, may exhibit such correlation. 
Both types of decays can have direct and mixing-induced CP violation. A third type of decays 
involves charged kaons and does not exhibit mixing-induced CP violation, but usually has a strong 
correlation between the flavor tag and the Dalitz plot coordinate (the kaon charge), because it 
consists of S-flavor eigenstates. The Dalitz plot PDFs are obtained from MC simulation and are 
described with the use of non-parametric functions. For neutral B background, the signal At 
resolution model is assumed. 

The Dalitz plot treatment of the continuum events is similar to the one used for charged-S 
background. The square Dalitz plot PDF for continuum background is obtained from on-rcsonance 
events selected in the rriES sidebands and corrected for feed-through from B decays. A large 
number of cross checks has been performed to ensure the high fidelity of the empirical shape 
parameterization. Analytical models have been found insufficient. The continuum At distribution 
is parameterized as the sum of three Gaussian distributions with common mean and three distinct 
widths that scale the At pcr-cvcnt error. This yields six shape parameters that arc determined by 
the fit. The model is motivated by the observation that the At average is independent of its error, 
and that the At RMS depends linearly on the At error. 

3.3.2 PARAMETERIZATION OF THE OTHER VARIABLES 

The mES distribution of TM signal events is parameterized by a bifurcated Crystal Ball func- 
tion [21], which is a combination of a one-sided Gaussian and a Crystal Ball function. The mean 
of this function is determined by the fit. A non-parametric function is used to describe the SCF 
signal component. 

The AE' distribution of TM events is parameterized by a double Gaussian function, where 
all five parameters depend linearly on rriQ. Misreconstructed events are parameterized by a broad 
single Gaussian function. 

Both mES and AE' PDFs are parameterized by non-parametric functions for all B-background 
classes. 

The TTiES and AE' PDFs for continuum events are parameterized with an Argus shape func- 
tion [23] and a second order polynomial, respectively, with parameters determined by the fit. 

We use non-parametric functions to empirically describe the distributions of the NN outputs 
found in the MC simulation for TM and SCF signal events, and for i?-background events. We 
distinguish tagging categories for TM signal events to account for diff'erences observed in the shapes. 

The continuum NN distribution is parameterized by a third order polynomial that is defined to 
be positive. The coefficients of the polynomial are determined by the fit. Continuum events exhibit 
a correlation between the Dalitz plot coordinate and the shape of the event that is exploited in the 
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NN. To correct for residual effects, we introduce a linear dependence of the polynomial coefficients 
on the distance of the Dalitz plot coordinate to the kinematic boundaries of the Dalitz plot. The 
parameters describing this dependence are determined by the fit. 

4 SYSTEMATIC STUDIES 

The contributions to the systematic error on the signal parameters are summarized in Table 3. 
Table 4 summarizes the correlation coefficient extracted from the systematic covariance matrix. 

The most important contribution to the systematic uncertainty stems from the signal modeling 
of the Dalitz plot dynamics. 

To estimate the contribution to -B° — > tt+tt^tt^ decay via other resonances and non-resonant 
decays, we have performed an independent analysis where we include these other decays in the fit 
model. For simplicity, wc assume a uniform Dalitz distribution for the non-resonance events and 
consider possible resonances including /o(980), /2(1270), and a low mass s-wave a. The fit does 
not find significant number of any of those decays. However, the inclusion of a low mass 7r"'"7r~ 
s-wave component degrades our ability to identify p°7r° events significantly. The systematic effects 
(contained in the "Dalitz plot model" field in Table 3) is estimated by observing the difference 
between the true values and Monte Carlo fit results, in which events are generated based on the 
new fit results and fit with the nominal setup where only p is taken into account. 

We vary the mass and width of the |o(770), p(1450) and p(1700) within ranges that exceed twice 
the errors found for these parameters in the fits to r and e^e~ data [11], and assign the observed 
differences in the measured U and / coefficients as systematic uncertainties {^^p^p' ^p" lineshape" in 
Table 3). Since some of the U and / coefficients exhibit significant dependence on the p(1450) and 
p(1700) contributions, we leave their amplitudes (phases and fractions) free to vary in the nominal 
fit. 

To validate the fitting tool, we perform fits on large MC samples with the measured proportions 
of signal, continuum and S-background events. No significant biases are observed in these fits. The 
statistical uncertainties on the fit parameters are taken as systematic uncertainties ("Fit bias" in 
Table 3). 

Another major source of systematic uncertainty is the i?-background model. The expected 
event yields from the background modes are varied according to the uncertainties in the measured 
or estimated branching fractions ("-/VBackground" in Table 3). Since 5-background modes may exhibit 
CP violation, the corresponding parameters are varied within appropriate uncertainty ranges ("B 
background CP" in Table 3). As is done for the signal PDFs, we vary the At resolution parameters 
and the flavor-tagging parameters within their uncertainties and assign the differences observed in 
the on-resonance data fit with respect to nominal fit as systematic errors. 

Other systematic effects are much less important to the measurements of U and / coefficients, 
and they are combined in "Others" field in Table 3. Details are given below. 

The parameters for the continuum events are determined by the fit. No additional systematic 
uncertainties are assigned to them. An exception to this is the Dalitz plot PDF: to estimate 
the systematic uncertainty from the twes sideband extrapolation, wc select large samples of off- 
resonance data by loosening the requirements on Ai? and the NN. We compare the distributions 
of m! and 9' between the mES sideband and the signal region. No significant differences are found. 
We assign as systematic error the effect seen when weighting the continuum Dalitz plot PDF by 
the ratio of both data sets. This effect is mostly statistical in origin. 
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The uncertainties associated with Am^i and r are estimated by varying these parameters within 
the uncertainties on the world average [20]. 

The systematic effects due to the signal PDFs comprise uncertainties in the PDF parameter- 
ization, the treatment of misreconstructed events, the tagging performance, and the modeling of 
the signal contributions. 

When the signal PDFs are determined from fits to a control sample of fully reconstructed B 
decays to exclusive final states with charm, the uncertainties are obtained by varying the parameters 
within the statistical uncertainties. In other cases, the dominant parameters have been left free to 
vary in the fit, and the differences observed in these fits are taken as systematic errors. 

The average fraction of misreconstructed signal events predicted by the MC simulation has been 
verified with fully reconstructed B Dp events [17]. No significant differences between data and 
the simulation were found. We vary /scF ^oi all tagging categories relatively by 25% to estimate 
the systematic uncertainty. 

Tagging efficiencies, dilutions and biases for signal events are varied within their experimental 
uncertainties. 

The systematic errors for the parameters that measure interference effects are dominated by 
the uncertainty in the signal model, mainly the tail description of the p resonance. For the other 
parameters, the uncertainty on the fit bias and the S-background contamination are important. 

5 FIT RESULTS 

The maximum-likelihood fit results in the B^ 7r^-7r~7r'^ event yield 1847 it 69, where the error is 
statistical only. For the U and / coefficients, the results are given together with their statistical and 
systematic errors in Table 5. The corresponding correlation matrix is given in Table 6. We have 
generated a sample of Monte Carlo experiments to determine the probability density distributions of 
the fit parameters. Within the statistical uncertainties of this sample we find Gaussian distributions 
of the distribution for the fitted U and / coefficients. This allows us to use the least-squares method 
to derive other quantities from these (Section 6). 

The signal is dominated by B^ — > p^ir^ decays. We observe an excess of /O^'tt^ events, which 
is in agreement with our previous upper limit [24], and the latest measurement from the Belle 
collaboration [25]. The result for the p(1450) amplitude is in agreement with the findings in r and 
e+e~ decays [11]. For the relative strong phase between the p(770) and the yo(1450) amplitudes we 
find (171 lb 23)° (statistical error only), which is compatible with the result from r and e^e~ data. 

Figure 3 shows distributions of AE', m-Es, the NN output, At/a{At), where a{At) is the per- 
event error on At, as well as the Dalitz plot variables m' and 6', which are enhanced in signal 
content by requirements on the signal-to-continuum likelihood ratios of the other discriminating 
variables. Figure 4 shows distribution of the minimum of three di-pion invariant masses, again 
enhanced in signal content. This plot shows clearly that p(770) dominates the signal component. 

As a validation of our treatment of the time dependence we allow r^o to vary in the fit. We 
find Tgo = (1.513 lb 0.066) ps, while the remaining free parameters are consistent with the nominal 
fit. To validate the SCF modeling, we leave the average SCF fractions per tagging category free to 
vary in the fit and find results that are consistent with the MC prediction. 
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6 INTERPRETATION OF THE RESULTS 



The U and / coefficients are related to tlie quasi-two-body parameters, defined in Ref. [17] as 
follows 

^^ = ^' ^~ = ^' ^^ = lJf^ ^~^lJf' ^'^- = u}+u} ' ^^^^ 

and where C = {C+ + C-)/2, AC = (C+ - C-)/2, S = {S+ + S-)/2, and AS = {S+ - S-)/2 
and ApTT is the time and flavor integrated asymmetry. In contrast to our previous analysis [17], the 
definitions of Eq. (26) explicitly account for the presence of interference effects, and are thus exact 
even for a p with finite width, as long as the U and / coefficients are obtained with a Dalitz plot 
analysis. This treatment leads to a dilution of the result and hence to slightly increased statistical 
uncertainties compared to neglecting the interference effects. 
For the CP-violation parameters, we obtain 

Apn = -0.142 ± 0.041 ± 0.015 , (27) 
C = 0.154 ± 0.090 ± 0.037 , (28) 
S = 0.01 ± 0.12 ± 0.028 , (29) 

where the first errors given are statistical and the second are the systematic uncertainties. For the 
other parameters in the quasi-two-body description of the B'^{B^) pir decay-time dependence, 
we measure 

AC = 0.377 ± 0.091 ± 0.021 , (30) 
AS = 0.06 ± 0.13 ± 0.029 . (31) 

The systematic errors are dominated by the uncertainty on the CP content of the B-related back- 
grounds. Other contributions are the signal description in the likelihood model (including the 
limit on non-resonant vr+vr^vr'^ events), and the fit bias uncertainty. The covariance matrix, 

including systematic, of the five quasi-two-body parameters is given in Table 7. 

One can transform the experimentally convenient, namely uncorrelated, direct CP-violation 
parameters C and Ap^ into the physically more intuitive quantities A'^, A'^, defined by 



^ k+~P-l ^ Ap^ + C + Ap^AC 
-^f^ |k+-|2 + 1 l + AC + ^p,C ' ^ ' 



A — I- _ 1^ 1 ■^(m C ApjiAC 



+ 1 1 - AC - Afy„C 



where = {q/p){A 7^4+) and k + = (g/p)(A+/yl ), so that .4^ (.4^^) involves only diagrams 
where the p (tt) meson is emitted by the W boson. We find 

At~ = 0.03 ± 0.07 ± 0.03 , (33) 



^p7V 

A-^ = -0.38l°;i|±0.07, (34) 



with a correlation coefficient of 0.62 between and Ap^. The confidence level contours including 
systematic errors are given in Fig. 5. The significance, including systematic uncertainties and 
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calculated by using a mininum method, for the observation of non-zero direct CP violation is at 
the 3.0(7 level. The evidence of direct CP violation is almost entirely from the p~'K^ decays. 

The measurement of the resonance interference terms allows us to determine the relative phase 

5+_=arg(^+M-) , (35) 

between the amplitudes of the decays —>■ p~Tr~^ and B^ — > p+Tr". Through the definitions (8)- 
(12), we can derive a constraint on S+- from the measured U and / coefficients by performing a 
least-squares minimization with the six complex amplitudes as free parameters. The constraint 
can be improved with the use of strong isospin symmetry. The amplitudes A'^ represent the sum 
of tree-level and penguin-type amplitudes, which have different CKM factors: the tree-level {T'^) 
_^ pK^K transition amplitude is proportional to VudV*f^, while the corresponding penguin-type 
amplitude {P'^) comes with VqdV*^^, where q = u,c,t. Here we denote by k the charge conjugate 
of K, where = 0. Using the unitarity of the CKM matrix one can reorganize the amplitudes and 
obtain [12] 

^« ^ j^K^-ia ^ pK ^ 

where the magnitudes of the CKM factors have been absorbed in the T"^, P", T'^ and P'^. The 
Eqs. (36) represent 13 unknowns of which two can be fixed due to an arbitrary global phase and the 
normalization condition = 1. Using strong isospin symmetry and neglecting isospin-breaking 
effects, one can identify P° = —{P^ + P~)/2, which reduces the number of unknowns to be 
determined by the fit to 9. This set of parameters provides the constraint on shown in the 
left plot of Fig. 6. We find for the solution that is favored by the fit 

(5+_ = (34 ±29)° , (37) 

where the errors include both statistical and systematic effects. There is a clear structure of multiple 

solutions which give comparable x^. Only a marginal constraint on 5-\ is obtained for C.L. < 0.05. 

Finally, following the same procedure, we can also derive a constraint on a from the measured U 
and I coefficients. The resulting C.L. function versus a is given in the right hand plot of Fig. 6. It 
includes systematic uncertainties. Ignoring the mirror solution at a + 180°, we find a G (75°, 152°) 
at 68% C.L. No constraint on a is achieved at two sigma and beyond. 

7 SUMMARY 

We have presented the preliminary measurement of CP-violating asymmetries in B^ vr+Tr^vr^ 
decays dominated by the p resonance. The results are obtained from a data sample of 347 million 
T(45) BB decays. We perform a time-dependent Dalitz plot analysis. From the measurement 
of the coefficients of 26 form factor bilinears we determine the three CP-violating and two CP- 
conserving quasi-two-body parameters, where we find a 3.0fj evidence of direct CP violation. Taking 
advantage of the interference between the p resonances in the Dalitz plot, we derive constraints on 
the relative strong phase between B^ decays to p+Tr" and p^tt"*", and on the angle a of the Unitarity 
Triangle. These measurements are consistent with the expectation from the CKM fit [26]. 
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Figure 1: Jacobian determinant (17) of the transformation (15) defining the SDP. Such pattern 
would be obtained in the SDP if events were uniformly distributed over the nominal Dalitz plot. 




5 10 15 20 25 30 0.2 0.4 0.6 0.8 1 



(GeVV) m' 

Figure 2: Nominal (left) and square (right) Dalitz plots for Monte-Carlo generated — > tt^tt^tt^ 
decays. Comparing the two Dalitz plots shows that the transformation (15) indeed homogenizes 
the distribution of events which are no more along the plot boundaries but rather cover a larger 
fraction of the physical region. The decays have been simulated without any detector effect and the 
amplitudes A" and have all been chosen equal to 1 in order to have destructive interferences 
at equal p masses. The main overlap regions between the charged and neutral p bands are indicated 
by the hatched areas. Dashed lines in both plots correspond to ^s+-fi = 1.5 GeV/c^: the central 
region of the Dalitz plot (defined by requiring that all 3 two-body invariant masses exceed this 
threshold) contains almost no signal event. 
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Class 


Mode 


BR [10-«] 


Expected number of events 





7->_L 1 n 


19.1 ±3.5 


52 ± 10 





(— (/?7r)+)7r" 


20.0 ±15.0 


32 ± 24 





7->_|_ n / -A U\ 4- 

B+ 0-1 p TT +)7r+ 


20.0 ± 15.0 


19 ± 14 


1 




8.7 ± 


1.0 


73 ± 8 


1 


B^ —>■ p 


4.3 ± 


0.6 


C 1 1 

6 ± 1 


2 


B^ — 7r^iVg(— ^ TT^TT ) 


8.3 ± 


0.4 


10 ± 1 


6 


D+ + 

B^ TT p^ 


10.8 ± 


1.4 


63 ± 8 


Q 
O 


R+ ^ Tr+ 1^0 1. ^0_0\ 

±) — > TT J^g{ — > TT TT ) 


3.7 ± 


0.2 


it) ± z 


A 


-D ? TT TT 


5.5 ± 


0.6 


14-1- 9 


A 

4 


— ^ A ^TT 


12.1 ± 


0.8 


_L 1 





D+ V ZV^r**) /I /I -T^\ + V f Tf^-^^^X-i- 


29.0 ± 


5.4 


OQ _1_ K 

6o ± 





B" TT A^(— >■ KgTr^) 


3.3 ± 


0.4 


2 ± 1 


IT 
i 


^ ^^>[long] 


25.2 ± 


3.7 




7 


B" — s- (oitt) 


39.7 ±3.7 


39 ± 4 


8 


B" — ^ iV^TT 


18.9 ± 


0.7 


12 ± 


9 


B IT K K^TT ) 


3.3 ± 


0.4 


20 ± 2 


9 


B" a'- ''(1430)7r Ktttt^ 


11.2 ± 


2.2 


212 ± 34 


iU 


Jd — > 'yJ\ {ovZ, i4oUJ(^ — > [J\ TT ) ) 


27.4 ± 


1.5 


1 /I -L 1 


1 n 


t> — ^ TT A ( — > A ' TT j 


1.3 ± 


0.5 


y ± 4 


1 n 


£j ^ 1] p 7j7r 


0.4 ± 


0.2 


o HZ z 


11 




9.9 ± 


1.6 


103 ± 17 


12 


[nonresj 


4.6 ± 


4.6 


38 ± 38 


13 


fiO ^ TT+TT-) 


5.8 ± 


0.5 


50 ± 4 


14 


^0 ^ Z)-(^ 7r-7r0)7r+ 


7.5 ± 


2.3 


599 ± 184 


15 


^0 ^ ^0(^ K+vr-)7r0 


11.0 ± 


3.2 


100 ± 29 


16 


^0 ^ ^0(^ vr+7r-)7rO 


0.4 ± 


0.1 


35 ± 9 


17 




2.6 ± 


0.5 


77 ± 15 


18 


B^ {neutral generic b ^ c decays} 






173 ± 15 


19 


S+ {charged generic b ^ c decays} 






396 ± 20 



Table 1: Summary of the S-background modes taken into account for the likelihood model. They 

have been grouped in twenty classes: charged charmless (six), neutral charmless (eight), exclusive 
neutral charmed (four) and inclusive neutral and charged charmed decays. Modes with at least 
two events expected after final selection have been included. 
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Variable 


TM Signal 


SCF Signal 


Continuum 


B-Background 


AE 


GG 


G 


P2 


NP 


rriES 


biCB 


NP 


Argus 


NP 


Neural Net 


NP 


NP 


P3 


NP 


Dalitz 


see text 


sec text 


NP 


NP 


At 


Ct Ct Ct 


Ct Ct Ct 


Ct Ct Ct 


Ct Ct Ct 



Table 2: Summary of PDF parameterizations where G=Gaussian, PX=X-order polynomial, 
NP=non-parametric, and biCB=bifurcated Crystal Ball. See Section 3.3.1 for a detailed description 
of the Dalitz plot parameterization for signal. 
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/n 


/_ 




—0 


—0 


/, 


+0 


Dalitz plot model 


0.010 


0.006 




0.110 


0.102 


0.020 


0.018 


p-,p' iP" linesliape 


0.003 


0.012 




0.240 


0.103 


0.009 


0.225 


Fit bias 


0.014 


0.023 




0.173 


0.375 


0.008 


0.186 


-^Background 


0.002 


0.008 




0.064 


0.085 


0.005 


0.026 


B background CP 


0.005 


0.009 




0.082 


0.061 


0.011 


0.039 


Others 


0.002 


0.006 




0.104 


0.052 


0.005 


0.051 


Sum 


0.017 


0.027 




0.315 


0.402 


0.023 


0.292 






rRe 
-'■+-0 


rim 










r7~,Re 


Dalitz plot model 


0.017 


0.007 


0.127 


0.082 


0.041 


0.144 


0.209 


PiP' iP" lineshape 


0.308 


0.138 


0.306 


0.012 


0.012 


0.086 


0.159 


Fit bias 


0.301 


0.036 


0.048 


0.088 


0.001 


0.050 


0.087 


-^Background 


0.049 


0.178 


0.176 


0.002 


0.009 


0.034 


0.045 


B background CP 


0.042 


0.044 


0.095 


0.015 


0.002 


0.020 


0.073 


Others 


0.073 


0.059 


0.095 


0.004 


0.002 


0.012 


0.050 


Sum 


0.431 


0.143 


0.335 


0.121 


0.043 


0.175 


0.276 





r/-+,Im 


TT+,Re 


uz 






7-7— , Re 


7-7-+,Im 


Dalitz plot model 


0.034 


0.024 


0.022 


0.030 


0.036 


0.258 


0.076 


p^p\p" lineshape 


0.222 


0.045 


0.010 


0.030 


0.050 


0.216 


0.089 


Fit bias 


0.034 


0.058 


0.004 


0.005 


0.053 


0.007 


0.004 


-^Background 


0.045 


0.017 


0.013 


0.005 


0.117 


0.103 


0.034 


B background CP 


0.032 


0.018 


0.032 


0.009 


0.063 


0.075 


0.019 


Others 


0.027 


0.015 


0.010 


0.002 


0.045 


0.020 


0.027 


Sum 


0.227 


0.078 


0.025 


0.043 


0.082 


0.337 


0.117 





rr+,Re 




^-3e 










Dalitz plot model 


0.045 


0.014 


0.250 




0.703 


0.227 


0.010 


pjp' ,p" lineshape 


0.140 


0.169 


0.200 




0.169 


0.159 


0.031 


Fit bias 


0.020 


0.007 


0.069 




0.012 


0.033 


0.027 


-^Background 


0.069 


0.137 


0.122 




0.024 


0.166 


0.014 


B background CP 


0.014 


0.042 


0.029 




0.025 


0.025 


0.034 


Others 


0.024 


0.032 


0.067 




0.024 


0.044 


0.009 


Sum 


0.148 


0.170 


0.328 




0.723 


0.279 


0.042 








C 




AC 




S 




Dalitz plot model 


0.008 


0.002 




0.013 




0.015 


0.024 


p^p\p" lineshape 


0.011 


0.021 




0.011 




0.016 


0.008 


Fit bias 


0.004 


0.012 




0.015 




0.028 


0.011 


-^Background 


0.002 


0.005 




0.011 




0.002 


0.011 


B background CP 


0.003 


0.029 




0.006 




0.017 


0.005 


Others 


0.001 


0.007 




0.005 




0.006 


0.007 


Sum 


0.015 


o.oi^ 




0.021 




0.028 


0.029 



Table 3: Summary of systematic uncertainties. 
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-0.25 


-0.01 
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0.45 
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0.48 
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-0.56 


-0.53 
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-0.89 


-0.54 


0.19 


0.31 


-0.28 


-0.68 


0.12 


-0.04 


0.06 


-0.38 


-0.64 


0.88 


0.84 


j^+,Rc 


0.56 


0.20 


0.00 


0.13 


0.18 


0.70 


0.28 


-0.30 


-0.29 


0.00 


0.50 


-0.55 


-0.37 




-0.04 


0.56 


-0.54 


-0.44 


0.44 


-0.16 


-0.53 


0.65 
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0.30 
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-0.23 


-0.18 
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0.49 
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0.11 


1.00 
















7r+,Im 


-0.55 


0.17 


0.40 


0.27 


-0.52 


-0.10 


0.17 


1.00 
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-0.20 


-0.59 


0.48 


-0.16 


-0.37 


0.42 


-0.54 


0.10 


1.00 














-0.48 


-0.70 


0.10 


-0.19 


-0.39 


-0.53 


-0.56 


0.21 


0.23 


1.00 












0.48 


-0.20 


-0.30 


-0.51 


0.65 


-0.08 


-0.55 


-0.69 


-0.13 


0.21 


1.00 








^+,Im 


-0.68 


0.21 


0.30 


0.48 


-0.62 


-0.24 


0.68 


0.68 


-0.28 


0.09 


-0.72 


1.00 








0.20 


-0.62 


-0.19 


-0.47 


0.18 


-0.27 


-0.91 


-0.28 


0.35 


0.60 


0.61 


-0.68 


1.00 






0.20 


0.48 


-0.02 


0.34 


0.56 


0.30 


0.19 


-0.29 


-0.36 


-0.33 


0.31 


-0.17 


-0.14 


1.00 



Table 4: Correlation matrix of systematic uncertainties for the U and / coefficients. Note that all 
elements above the diagonal are omitted for readability. 
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Parameter Description Result 

Coefficient of |/op 0.237 ± 0.053 ± 0.043 
U+ Coefficient of 1.33 ± 0.11 ± 0.04 

Uq Coefficient of |/op cos(AmdAt) -0.055 ± 0.098 ± 0.13 

UZ Coefficient of cos(AmrfAt) -0.30 ± 0.15 ± 0.03 

[/+ Coefficient of cos(AmrfAt) 0.53 ± 0.15 ± 0.04 

lo Coefficient of |/op sin(AmdAt) -0.028 ± 0.058 ± 0.02 

/_ Coefficient of sin(AmdAi) -0.03 ± 0.10 ± 0.03 

/+ Coefficient of sin(AmdAt) 0.039 ± 0.097 ± 0.02 

CZ+i^"" Coefficient of Im[/+/*] 0.62 ± 0.54 ± 0.72 

Coefficient of Re[/+/*] 0.38 ± 0.55 ± 0.28 

Coefficient of Im[/+/*] cos(AmrfAt) 0.13 ± 0.94 ± 0.17 

CZ+i^'' Coefficient of Re[/+/*] cos(AmrfAt) 2.14 ± 0.91 ± 0.33 

4™. Coefficient of Im[/+/* ] sin(AmdAt) -1.9 ±1.1 ±0.1 

Ifl Coefficient of Re[/+/*] sin(AmrfAt) -0.1 ±1.9 ±0.3 

U:^^'^ Coefficient of Im[/+/o*] 0.03 ± 0.42 ± 0.12 

U:!:^ Coefficient of Re[/+/o*] -0.75 ± 0.40 ± 0.15 

U~^'^ Coefficient of Im[/+/*] cos(Am<iAt) -0.93 ± 0.68 ± 0.08 

U'^^"" Coefficient of Rc[/+/o*] cos(AmdAt) -0.47 ± 0.80 ± 0.3 

I^^ Coefficient of Im[/+/o*] sin(AmdAt) -0.1 ±1.1 ±0.3 

7^3 Coefficient of Rc[/+/*] sin(AmdAt) 0.2 ±1.1 ±0.4 

[7+^^"^ Coefficient of Im[/_/o*] -0.03 ± 0.40 ± 0.23 

[7+^^*" Coefficient of Re[/_/o*] -0.52 ± 0.32 ± 0.08 

UZ^'^ Coefficient of Im[/_/o*] cos(AmrfAt) 0.24 ±0.61 ±0.2 

UZ^ Coefficient of Re[/_/o*] cos(AmdAt) -0.42 ± 0.73 ± 0.28 

I^Z^ Coefficient of Im[/_/o*]sin(Amrf At) 0.7 ±1.0 ±0.3 

I^^ Coefficient of Re[/_/o*] sin(AmdAt) 0.92 ± 0.91 ±0.4 

Table 5: Fit results for the U and I coefficients. The errors given are statistical (first) and systematic 
(second). The free normalization parameter is fixed to 1. 
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Table 6: Correlation matrix of statistical uncertainties for the U and / coefficients. Note tliat all 
elements above the diagonal are omitted for readability. 
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Figure 3: Distributions of (top to bottom, left to right) /S.E\ tues, NN output, At/a{At), m' and 
9' for samples enhanced in it'^tt~7t^ signal. The dots with error bars give the on-resonance 

data. The solid histogram shows the projection of the fit result. The dark, medium and light 
shaded areas represent respectively the contribution from continuum events, the sum of continuum 
events and the i?-background expectation, and the sum of these and the misreconstructed signal 
events. 
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Figure 4: Distribution of minimum of the three di-pion invariant masses, for samples enhanced in 
— > 7r'*'7r~7r'^ signal. The dots with error bars give the on-resonance data. The solid histogram 
shows the projection of the fit result. The dark, medium and light shaded areas represent respec- 
tively the contribution from continuum events, the sum of continuum events and the ^-background 
expectation, and the sum of these and the misreconstructed signal events. 
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Table 7: Covariance matirx multipled by 1000 of the quasi-two-body parameters. 
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Figure 5: Confidence level contours for the direct CP asymmetries A'^^ versus Ap^^. The shaded 
areas represent la, 2a and 3a contours, respectively. 
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Figure 6: Confidence level functions for 5^ (left) and a (right). Indicated by the dashed horizontal 

lines are the confidence level (C.L.) values corresponding to la and 2a, respectively. 
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